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We have investigated theoretically the interaction between individual quantum dot with broken 
inversion symmetry and electromagnetic field of a single-mode quantum microcavity. It is shown 
that in the strong coupling regime the system demonstrates nonlinear optical properties and can 
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I. INTRODUCTION 

Quantum microcavity is unique laboratory for stud- 
ies of the strong light-matter coupling. Being first ob- 
served two decades ago [l| , the strong coupling regime is 
now routinely achieved in different kinds of microcavities 
[^. From the fundamental viewpoint, it is interesting 
as a basis to investigate various collective phenomena in 
condensed matter systems such as Bose-Einstein conden- 
sation (BEG) 3 and superfluidity Q. From the view- 
point of applications, it opens a way to the realization 
of optoelectonic devices of the new generation [5| : room- 
temperature polariton lasers [6[, polarization- controlled 
optical gates [3], and others. 

Several applications of the strong coupling regime were 
also proposed for the quantum information processing 
[8l-[l0l|. In this case one should be able to tune the num- 
ber of emitted photons in controllable way. This is hard 
to achieve in planar microcavities where the number of 
elementary excitations is macroscopically large, but is 
possible in microcavities containing single quantum dots 
(QDs) where QD exciton can be coupled to a confined 
electromagnetic mode provided by a micropillar (etched 
planar cavity) [il[ , a defect of the photonic crystal [l^] or 
a whispering gallery mode [l^- That is why the strong 
coupled systems based on QDs have attracted particular 
attention recently. 

Depending on the size of the QD. its elementary exci- 
tations can behave as fermions (small QDs whose size is 
comparable with exciton radius in the bulk material) [l4| , 
bosons (large QDs whose size is much larger than exciton 
radius in the bulk material) |l5| or particles with interme- 
diate statistics (medium size QDs) [lal- In the current 
paper, we consider a small single QD in a microcavity 
that corresponds to the case of fermions. For symmet- 
ric QDs such system can be described by the well-known 
Jaynes-Cummings Hamiltonian [17{ which predicts trans- 
formation of the Rabi doublet to the Mollow triplet in 
the emission spectrum as intensity of the external pump 
growth for both coherent [T^l and incoherent excitation 
schemes [13- On the other hand, the incorporation of 
the asymmetry into the quantum system can radically 
change its emission pattern and lead to the opening of 
optical transitions which were forbidden in the symmet- 
ric case. Particularly, the breaking of inversion symmetry 
opens optical transitions at the Rabi frequency at QDs 
placed in strong external laser field [2y|. The similar ef- 



fect occurs for asymmetric quantum wells placed inside 
a planar microcavity J2ll |22| . In current manuscript we 
consider modification of the emission spectrum of asym- 
metric QDs inside a single-mode microcavity using fully 
quantum approach. 

The work is organized as follows. In Section II we 
describe the formalism and introduce the model Hamil- 
tonian. In Section III we obtain analytical solutions for 
the important particular cases. In Section IV we dis- 
cuss the incorporation of pump and decay terms into the 
Hamiltonian and present the numerical calculations of 
the emission spectrum. Section V contains discussions 
and conclusions. 



II. MODEL 

We model QD as a two-level quantum system with 
the ground state \g) and excited state |e) with energies 
Eg and Ee, respectively. QD is placed inside a cavity 
and interacts resonantly with confined electromagnetic 
mode of the frequency Wc. Since electromagnetic field can 
transfer an electron in the QD from the valence band into 
the conduction band, the ground state \g) corresponds to 
the absence of free carriers while the first excited state |e) 
is the state with an electron in the conduction band and a 
hole in the valence band. Therefore, the energy difference 
A ~ Eg ~ Ee is approximately equal to the band gap 
of the QD minus excitonic correction accounting for the 
Goulomb attraction between the electron and the hole. 

The full Hamiltonian of the system can be represented 
as a sum of three parts. 



where 



H — He + 'Huj + Hint 






(1) 



(2) 



is the Hamiltonian of the single QD, and (Jx,y,z are the 
Pauli matrices acting in the space of |e) and \g) states. 
The Hamiltonian of the free electromagnetic field reads 



H^ 



fujjf.a* a, 



(3) 



where a, a' are the annihilation and creation operators 
for cavity photons, respectively. The Hamiltonian "Hint 
describes interaction of the QD with the electromagnetic 



field and can be constructed as following. The inter- 
action of a classical dipole d with a classical external 
electric field E is given by the expression T-Liat ~ —Ed. 
Within the quantum-field approach, we have to replace 
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FIG. 1: Two-level quantum dot with the bandgap A in a 
single- mode microcavity with the frequency u)c- 

the classical quantities d and E with the corresponding 
operators 
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where e is the polarization vector of the cavity mode, eo 



is the vacuum dielectric permittivity, d^ 



{i\d\j) are 



dipole matrix elements of the QD, I is the unity matrix, 
(T^ = {ax ± iay)/2. In the case of symmetric QD, the 
matrix elements dee and Agg are zero. Due to the break- 
ing of inversion symmetry in asymmetric QDs, the dipole 
matrix elements appear to be nonequivalent, dee 7^ ^gg- 
This leads to new physical effects discussed hereafter. 
The interaction Hamiltonian can be written as 



Hint = -dE 

= fffl (a + at) (ct+ + a-) + gs {a + a'') {n, + I) 

« gB^ {aa+ + aV") +gs{a + a)) [a, + /) . (5) 

The coupling parameters, gji = ^ (dcg • e)-y/?ia;c/2eoT^ 
and gs = ^ (dee • e) ^/nw^jAeoV, describe the excitation- 
exchange interaction and the symmetry-dependent inter- 
action, respectively. For definiteness, we assume them to 
be real. The parameter V in the expressions above is the 
quantization volume and can be estimated as ^ « (A/2)'^, 
where A = c/2tiuJc is the characteristic wavelength cor- 
responding to the cavity mode. We also put dgg = 0, 
which can be justified for nonferroelectric QDs. Indeed, 
asymmetry of such QDs is provided by peculiar shape 
and/or an external electric field. Due to it, the matrix 
element dgg is proportional to the size of the elementary 
cell of the crystal lattice, while the matrix element dee is 
proportional to the size of the QD. As a result, in realis- 
tic QDs one has dee ^ dgg ■ To pass from the second line 
to the third line in Eq. (jS]), the rotating- wave approx- 
imation [23J, |2J] was applied and we dropped the anti- 
resonant terms proportional to aV^ and aa^ . 



The full Hamiltonian of the system reads 



n 



f. t ^ 



+gR (aCT+ + aV ) + gs {a + a^) (cr^ + /) . 



(6) 



In the case of a symmetric QD, the coupling parameter 
gs equals zero. Then Eq. ([5]) reduces to the Hamiltonian 
of the fully solvable Jaynes-Cummings model. Its eigen- 
states correspond to electronic excitations dressed by the 
cavity photons and can be expressed as 



|^±(°))=A±|5,n)+i3±|e,n-l), 



where 



-±(0) 



-ii-^i 



Bt = 



hbJcin - 1) - A/2 



[ef}'' 



hwcin - 1) - A/2]2 + 



9%^ 



9R\/n 



-±(0) 



tvjjc{n-l) - A/2] 2 + g|r 



(7) 



(8) 



(9) 



and the composite electron-photon states |5,ri) = 
1.9) <2) I") and |e, n) = |e) (g) |n) describe both the QD 
state (the ground state g or the excited state e) and the 
field state with n cavity photons. It should be noted 
that the Jaynes-Cummings Hamiltonian commutes with 
the excitation number operator 



N = a\L + {a,+I)/2, 



(10) 



whose eigenstates correspond to the conserved number of 
total electron-photon excitations in the system counted 
as number of the excitations in QD (one for the state 
\g), zero for the state |e)) plus number of the photons in 
the cavity mode. The eigenenergies corresponding to the 
states d?]) are given by 



.±(0) 



hWr 



^)-;^^^ 



9Rn- 



ill) 



In order to obtain the emission spectrum of the system, 
we need to analyze optical transitions between the eigen- 
states ([7]). An emitted photon goes outside the system, 
and a pumped photon appears in the system. Thus, 
there is an exchange of photons between the coupled QD- 
microcavity system and some external reservoir. There- 
fore, we can introduce the Hamiltonian of the exchange 
of photons between cavity and outside world. 



7ie:E ^h{T*ar^ +Ta)r), 



(12) 



where r ' , r are creation and annihilation operators for 
the external photons, F is the system-reservoir coupling 
constant. 

The probabilities (intensities) of transitions with emis- 
sion of photon from the cavity are proportional to the 
corresponding matrix elements. 



hf-\{i^fAR\'HeAi^.,OR)? 



(13) 



where the symbols ipf , ^i denote the final and initial 
eigenstates of the Hamiltonian (O, Or and 1r describe 



zero- and one-photon states of the reservoir. Substituting 
Eq. ((HI) into Eq. dH]), one gets 



^'^'if^r, 









(14) 



where the states \4'i.f) are defined by Eq. ([7|), and integers 
ni,nf are initial and final numbers of electron-photon 
excitations in the system, defined earlier (see Eq. ((TO| V 

The Kroneckcr delta in Eq. (fT4| means that only tran- 
sitions changing number of excitations by one are allowed 
in a system described by the Jaynes-Cummings Hamil- 
tonian. It follows from Eq. (|lip that in the case of reso- 
nance {hiOc = A) the optical spectrum contains peaks at 
the energies fluJc±gR{Vn + liv^); where n = 0, 1,2, .... 
The accounting of the peaks broadening leads to the 
emission spectrum in the form of the Rabi doublet (for 
n = 0), and the Mollow triplet (for n > 1) [ll]. In the in- 
termediate regime more complicated multiplet structure 
can be observed |14| . 

It should be noted that in symmetric QDs {gg = 0) 
transitions at the Rabi frequency ^Ir = 2gii/h are for- 
bidden. Indeed, the transitions would occur between 
electron-photon states with the same number of excita- 
tions, 71/ = Hi, that is not allowed by Eq. (|T4| . However, 
these transitions become possible for asymmetric QDs, 
i.e. when the full Hamiltonian ^ contains the term 
gs (a -I- a'^) {a^ + 1) [20] • Since for realistic microcavi- 
ties the Rabi frequency ^r lies in the terahertz (THz) 
range, such transitions form the physical basis for using 
the considered system as a tunable source of THz radi- 
ation. We will consider these THz transitions in more 
details in the following sections. 



III. ANALYTICAL SOLUTIONS 

First of all, let us consider analytically the case of 
weak asymmetry. Then the Hamiltonian (|6]) can be rep- 
resented as a sum of two parts 



n^Ujc + v, 



where 



T^.JC — ^ca a + -z-^z + gniaa^ + a' a ) 
is the Jaynes-Cummings Hamiltonian, and 

V :^gs{a + a)){a,+I) 



(15) 



(16) 



(17) 



is the term arising from the asymmetry of the QD. Con- 
sidering the term (J17p as a small perturbation and us- 
ing the standard first-order perturbation theory, the cor- 
rected wavcfunctions of the system can be written as 



iV'„^) = i^„^^°^) + EE 



K 



±a 



--±£ 



±(0) 



1(0) 



IC^°^ (18) 



It is easy to see that these matrix elements are different 
from zero only if tti = n ± 1, and thus states with n,n—l 
and 71 + 1 excitations become mixed by the asymmetry 
of QD. This mixture leads to the opening of the opti- 
cal transitions |V'n ) ~^ IV'n) at the frequencies U,f(y/n. 
As well, transitions at double frequencies of the cavity 
2ujc ± Vlji{y/n -f 1 ± y/n)/2 become opened. The allowed 
transitions are shown schematically in Fig. (2). 
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FIG. 2: Optical transitions in the asymmetrical QD coupled 
with the cavity mode. 

Since Rabi frequency Vtji lies typically in the THz 
range, asymmetric QD-cavity system can be used as a 
nonlinear THz emitter. Using Eqs. p^ and ([T5)) . we 
obtain the intensity of THz transitions as follows: 



+ \|2 



where 



hf-\{i^-\a\i^t)\ 



{i^n\a\^t) 



V-+ 

n—l,n 
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\/"-B„ B„+i) 



with coefficients A^ and Bf^ given by expressions ([5]) and 
([9]), respectively. 

One can also consider analytically another physical sit- 
uation: the asymmetry is no longer assumed to be a weak 
perturbation, but photon occupation numbers are sup- 
posed to be large, 7i 3> 1 . This limiting case corresponds 
to the classical electromagnetic field in the cavity. Let us 
represent the full Hamiltonian ([5]) as a sum of the "diag- 
onal" part 



where |V'n ) are the unperturbed eigcnfunctions ((T]), 
En are the corresponding eigenenergies (|11|) . and 



±a 



perturbation (fT7)) . 



/,"(o) 



("0" |^|'0m ) are the matrix elements of the 



iid = hhJcO^a + — 0-^ + gs{a + a^) {(^z + I) (19) 



and the "off-diagonal" part 



Tiod = gn (aa^ + aV ) 



(20) 



The Hamiltonian ([T9|) does not commute with the excita- 
tion number operator pU)) , but commutes with the Pauli 
matrix a^ ■ This means that eigenstates of the Hamilto- 
nian ([T^ can be represented as 



A.— 

oo 



(21) 
(22) 



fe=0 



After substituting the expressions ([21]) . p2|) into the 
Schrodinger equation "Hdj-i/'j'j''^) = £^''^l'0n'^) with the 
Hamiltonian ([T^ . we obtain the system of algebraic 
equations for coefficients C^jf : 



[fi^.fc - A/2 - £fj C^, = 



(23) 



for fc = 0,1,2,.., 



\hw,k + A/2 - <] Q, + 2g5Vfc+TQ,fc+i = (24) 
for fc = 0, 1, 

\hw,k + A/2 - el\ Q, 
+2.gs (\/^^TTCf,,,+i + VfcC;^,fc_i) = (25) 

for fc = 2,3,4,... 

The solutions of Eqs. ([25)1 are evident: 

A 



4 = fiwcn- — , 



(26) 



As for Eqs. (|25p . in the limiting case /c ;:^ 1 they are 
similar to the well-known recurrent expression for the 
Bessel function of the first kind: 



2mJ„i{x) = xJ„i-i{x) +xJ„i+i{x), 



(27) 



where m is an integer and x is the argument of the Bessel 
function of the first kind, J™ (x) . Therefore, we can write 
the solutions of Eqs. (PSJ) for fc ^ 1 as 



Then Xn = (dee • e)En/huJc, where En = y^2nhiJ^J7oV 
is the classical amplitude of the field. Therefore, in the 
case of the most relevant physical situation with a;„ ^ 1, 
the eigenfunctions ([5(7)) can be estimated as l?/)^) « \e,n) 
since Jfc_„(0) = S^^n- This allows to seek eigenfunctions 
of the full Hamiltonian 7i = Tid + 'Hod as a superposi- 
tion (linear combination) of the functions ([25)) and (|30l) . 
Substituting this superposition ifj^ into the Schrodinger 
equation H|^^) = e„ Jt/i^), wc can find the energy spec- 
trum e^ of the coupled electron-photon system. Par- 
ticularly, for the resonance case (A = huic) this linear 
combination is 



1^^ 



^i]rn)±\rn+,)), 



and the corresponding energies are 



hjJr 



■ee± 



Ml', 



(31) 



(32) 



where Vt'j^ = (deg • e)En/h is the Rabi frequency for the 
classically strong electromagnetic field. It follows from 
the expressions ([29)) - ([3T)) that the states with all pos- 
sible numbers of excitations become intermixed. This 
means that all transitions IV'n ) ~^ iV'm) ^-re allowed and 
emission spectrum contains frequencies {n — m)hwc^^']i- 
However, the intensity of the transitions decreases with 
increasing [n — m) because of decreasing the Bessel func- 
tions Jk{x) with increasing k [2^. Therefore, most in- 
tensive ones correspond to those depicted at Fig. [5] This 
agrees with the results obtained above in the frameworks 
of perturbation theory for the weak asymmetry case. 

It should be noted that there is no analytical solution 
for arbitrary photon occupation numbers of the cavity 
mode and arbitrary asymmetry strength. To study phys- 
ically relevant situations outside the obtained analytical 
solutions as well as to calculate the shape of the emis- 
sion spectrum, we need to apply the numerical approach 
discussed in the next section. 



IV. NUMERICAL APPROACH 



£^ = hujcn+ — , 

^nk = Jk-n[Xk) , 



(28) 



where Xk = —Ags\/k/fujJc. It should be noted that the 
solutions ([28)) satisfy Eqs. ([24 )) - ([25)) for smaU integers 
fc '--^ 1 as well, since 

lim Jn{x) — 0. 

n— »iboo 

As a result, for large photon occupation numbers n the 
eigenfunctions ([?T)) -([2 ^ take the form 



\rn 



\9,n) 



y^ Jfc-n(Xfc)|e, fc). 



(29) 
(30) 



fc=o 



Let us make n and V tend to infinity while keeping 
n/V constant. This limiting case corresponds to the con- 
ventional model of an intense laser-generated field |25j . 



The discussion in the previous section was dedicated 
to the analytical calculation of the energy spectrum of 
the system. In two limiting cases corresponding to weak 
asymmetry and large photon occupation numbers we 
were able to find the emission frequencies. However, even 
in these cases our treatment did not allow the calculation 
of the shape of the emission spectrum as function of the 
intensity of the external pump. In this section we cal- 
culate it numerically, using the approach based on the 
master equation techniques. 

Let us assume that one has a QD embedded in a micro- 
cavity and switches on incoherent pumping of the pho- 
tonic mode. After some time an equilibrium is estab- 
lished and steady state (SS) is reached. It means that the 
increase of photon number provided by external pumping 
is balanced by escape of photons from the cavity. In this 
regime, one can measure the emission spectrum, i.e. the 
intensity of the fiux of the photons going out from the 
cavity as a function of their frequency. 

In quantum optics, the standard way of consideration 
of the processes involving external pumping and decay is 



based on using of the master equation for the full density 
matrix of the system p (see, e.g., Ref . j2a| ) which can be 
represented in the following form: 



1 



dtP^ -r[U]p\ + Cp. 



(33) 



The first term on the right-hand side of the equation, H, 
stands for the Hamiltonian describing coherent processes 
in the system, and the symbol C denotes the Lindblad 
superoperator accounting for pump and decay. In the 
case we consider, T-L is given by the expression ([6]), while 
the Lindblad term reads: 

Cp = P [apa^ + a} pa — a* ap — paa^^ 
+ — ^ (2apa^ — pa} a — a^ap) 



2f i2apa 



Here P is the intensity of the incoherent pump of the cav- 
ity mode, "fph and jqd are broadenings of photonic and 
cxcitonic modes, respectively (the latter is taken zero in 
all calculations below) . Equation ([55)1 represents a set of 
linear ordinary differential equations for matrix elements 
of the density matrix p. In our numerical analysis we will 
use the basis 



\g,n), |e,n) 



(35) 



that gives us a following system of equations, which can 
be written briefly as 



dtpt' = MtJ^l^pli 



(36) 



.+ 



- pa^ a — a^ ap^ . (34) or, in the explicit form, as 






+ 9B. [VtTT6aeP-ll,j + VtSagptuj - Vl^bgptj-l " VJ^^beptj + l) 

+ gS [^/i + IfJaeP^+l J + ViSaePi-lJ - VjSbepfj-i - \/JTT5beptj + lj ' 2 

+ I (2V(* + l)(.? + l)Pm,,+i + 2 V^jpfi,,-! - (2* + 2j + 2)p^'^ 

+ ^-f (2v/(7TT)(7 + T)p.^l,, + l - i^+J)f^!) + ^ {'^ptM^b, - pf:,5ae - P^he) 



(37) 



Here the superscripts a, 6, c, d are either g or e (ground 
or excited state of the QD) and subscripts i,j,k,l cor- 
respond to the number of the photons in a cavity. In 
principle, they can take values from zero to infinity, but 
for numerical analysis truncation of the matrix is needed. 
The stronger is the pump, the more states should be 
taken into account. The natural way to control the ac- 
curacy of the truncation is to check the conservation of 
the trace of the truncated density matrix. 



Numerical solution of the system l\'3Q\i allows to find 
the full density matrix of the system in stationary state, 
pff , which allows to determine the probabilities of the 
occupancies of different quantum states of the coupled 
QD- cavity system as functions of the intensity of the 
external pump. As well, it allows to find the shape of the 
emission spectrum of the system. To pursue this latter 
task, we will use the two approaches. Let us start with 
the relatively simple model based on the modified Fermi 
golden rule, considering the isolated QD- cavity system. 
Its eigenstates can be found by diagonalization of the 
Hamiltonian ^ . This procedure was performed analyti- 
cally in the previous section for the cases of weak asym- 
metry and large photonic occupation numbers. However, 
in general case a numerical analysis is needed. Let the 
system is in a pure state corresponding to one of its eigen- 
states. Then, according to the Fermi golden rule, we can 



estimate the emission spectrum as 



(38) 



Tph 



f (£/ -e^- hwY + Yph 



|ai/Pp(^) 



llh 



{sf -e^- hujf + 7^^ ' 



where the symbols ipf , ij^i denote the final and initial 
eigenstates of the Hamiltonian ([5]), 0^, and 1^ describe 
zero- and one-photon states of the reservoir, T-Lex is the 
Hamiltonian of the coupling between the cavity and reser- 
voir (see EqfT2)). p(w) is the density of states in the reser- 
voir, aif — {tpi\a\^f). The Lorcntzian factor accounts for 
the broadening of the state, 7^/1, provided by finite life- 
time of cavity photons. 

If the system is in the mixed state, the spectrum ([38)) 
can be estimated as a sum over all possible initial states 
taken with corresponding probabilities Pi, 



S{u) 



27r 

T 






''p{^) 



ilh 



(£/ 



hwY 



ilh 



(39) 



Parameters Pi entering into the above expression are 
nothing but diagonal matrix elements of the density ma- 
trix in the stationary state written in the basis of eigen- 
states of the Hamiltonian ^, Pi — pfi^ . The matrix 
p can be found by an unitary transformation of the 



stationary density matrix in the basis (|e,ri,), \g,n)) ob- 
tained from solution of Eqs. ([551 ) 

The described phenomenological approach for calcu- 
lating the emission spectrum has one substantial draw- 
back. Namely, we assume the peaks in the spectrum to be 
Lorentzians which is not always guaranteed [2J, [23, l23 ■ 
In g eneral, according to the Wiener-Khintchine theorem 
[29j , the spectrum of the emission from the system can be 
calculated as a Fourier transform of two-time correlator. 



S{uj)r^ lim7^ / {a^{t + T)a{T))t 



'dr. 



(40) 



Generally, the calculation of two-time correlators is a 
complicated task which cannot be solved exactly. How- 
ever, it is well-known from literature [23, [2^ that the 
using of certain general assumptions about the behavior 
of the system allows to reduce the calculation of two-time 
correlators to calculation of one-time correlators within 
the framework of the the quantum regression theorem 
(QRT) [33 ■ This means that the spectrum can be cal- 
culated straightforwardly from the density matrix of the 
system in stationary state. For the system we consider 
in the present paper, the approach based in QRT results 
in 

5(c.) = -7^ Y. m+^f^^r']^]%iptl''^)^^r..a,.' 

7T ^ — ' ■'' 

i,j,k,l;a,b 

(41) 
where M is the matrix defined in Eqs. ([36[) . ([37[) . / is 
the unity matrix. ly should be noted that the expression 
above is valid for any choice of the basis. 

The shape of the emission spectrum of the system cal- 
culated using the Fermi golden rule and QRT is analyzed 
in the following section. 
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FIG. 3: Emission spectrum of the asymmetric QD- cavity 
system calculated using the Fermi golden rule. In optical 
diapason close to the eigenfrequency of the cavity standard 
quadruplet structure is revealed. In the range of frequencies 
close to the Rabi frequency additional set of emission peaks 
appears. They are provided by asymmetry of the QD and are 
absent in a symmetric case. The parameters of the calculation 
are: 7ph = 7150 =0.1 meV; P ;^ 0.15 meV; huc = A = 1 eV; 
Sfi w Ss « 1 meV. 
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V. DISCUSSION 

The emission spectrum of the system, calculated us- 
ing the phenomenological approach based on the Fermi 
golden rule, is presented in Fig. [3 

One sees that in the region of the frequencies close to 
the eigen frequency of the cavity (coinciding with the fre- 
quency of optical transition in QD, hbjc = A) the spec- 
trum reveals a quadruplet pattern. The appearance of 
the additional multiplets is not visible in the linear scale, 
but becomes apparent in logarithmic scale as it is seen 
in Fig. [H These results are in good qualitative agree- 
ment with those obtained earlier for the symmetric QDs 
in strong coupling regime [IJ]. This is not surprising, 
since main difference in the emission from symmetric and 
asymmetric QDs appears at the regions around the Rabi 
frequency and the double frequency 2a;c. 

The insets in Figs. I3I4I show the emission pattern at 
THz range about i^u. In full agreement with results of 
the Section II, one sees the appearance of the peaks in the 
emission at frequencies a; = flji^yn with n = 1,2,3,.... 
It should be noted that one should expect very low in- 
tensities of THz emission as compare to the emission in 
the optical diapason, since density of states of the photon 
reservoir scales as ut^. However, the situation can be im- 
proved if the coupled QD-cavity system is placed inside 
a bigger cavity tuned at the THz range. In this case the 



FIG. 4: Emission spectrum of the asymmetric QD- cavity 
system calculated using Fermi golden rule in logarithmic scale. 
Multiplet structure of the emission around eigenfrequency of 
the cavity becomes more visible then in the linear scale. In 
the range of the frequencies close to Rabi frequency additional 
set of emission peaks appears. Parameters of the system are 
the same as for Fig. [3| 



density of states has a sharp peak around the eigenfre- 
quency of THz cavity and the rate of spontaneous emis- 
sion is dramatically increased due to the Purcell effect 
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Current state of technology allows the increasing of 
the emission rates in the THz regime by a factor of 100 
in high-quality cavities [S^] ■ 

The results of the calculation of the spectrum based 
on using the quantum regression theorem are presented 
in Fig. [5j The part of the spectrum corresponding to the 
optical diapason for the frequencies about w « Wc is in 
good agreement with results obtained by using the Fermi 
golden rule (Figs. [31 [1|). The differences are that the side 
peaks in the QRT plots have lower intensities, and the 
shape of the spectrum is more smooth, so that multiplet 
pattern becomes invisible. 

The changes in the THz part are more dramatic. In- 
stead of series of the narrow peaks shown at the inset of 
Fig. [S]one sees a formation of a single broad strongly as- 
symmetric peak centered at a frequency ui « flR. How- 
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FIG. 5: Emission spectrum of the asymmetric QD- cavity 
system calculated using QRT. In optical diapason close to the 
eigenfrequancy of the cavity standard quadruplet structure 
is revealed. In the range of the frequencies close to Rabi 
frequency a single broad asymmetric emission peak appears. 
The parameters of the calculation are the same as in Fig. [3] 



ever, qualitative result remains the same: new optical 
transitions at THz range are opened by the asymmetry 
of the QD. 

The important question is statistical properties of the 
emitted THz light. As emission spectrum we obtain is 
quite large, and efficiency of THz emission is normally 
suppressed as compare to the emission at optical frequen- 
cies, one can expect that emitted radiation will have ther- 
mal statistics with second order coherence g'^' w 2. On 
the other hand, placing of the sample into high quality 
THz cavity can lead to selection of more narrow region 



of the frequencies of the emission. In this case one can 
expect the possibility to achieve THz lasing regime with 
5^^) « 1. The detailed consideration of this situation, 
however, lies beyond the scopes of the present paper. 
As to detection of the emitted THz radiation, it can be 
achieved by standard THz detectors (see, e.g., the review 



VI. CONCLUSIONS 

In conclusion, we considered an asymmetric two-level 
quantum system corresponding to small asymmetric QD 
interacting with an electric field of a single-mode micro- 
cavity. We found analytical solutions for the eigenener- 
gies of the system for the cases of weak asymmetry and 
high photon occupation numbers. As well, we developed 
numerical approach for the calculation of the emission 
spectrum under incoherent pumping. It is shown that in 
the regime of strong pump a new set of peaks in the emis- 
sion appears in the regions close to the Rabi frequency 
and double transition frequency. This allows to use the 
proposed system as nonlinear optical element and tun- 
able source of the THz radiation. 
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